Exact Solitonic Solutions of the Gross-Pitaevskii Equation with a Linear Potential 
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We derive classes of exact solitonic solutions of the time-dependent Gross-Pitaevskii equation 
with repulsive and attractive interatomic interactions. The solutions correspond to a string of 
bright solitons with phase difference between adjacent solitons equal to n. While the relative phase, 
width, and distance between adjacent solitons turn out to be a constant of the motion, the center 
of mass of the string moves with a constant acceleration arising from the inhomogeneouty of the 
background. 

PACS numbers: 



Introduction — The experimental realization of dark 
solitons 1 , 0, 01 , bright solitons [3, [H , and recently gap 
solitons 6] in Bose-Einstein condensates has stimulated 
intense interest in their properties particularly their for- 
mation and propagation 0, H, S 15, EH, 17 1. Due to 
the nonlinearity arising from the interatomic interactions 
and due to the presence of a confining potential in the 
Gross-Pitaeviskii equation that describes the evolution 
of the solitons, these studies were performed either by 
solving the corresponding Gross-Pitaevskii equation nu- 
merically or by using perturbative methods. Interest- 
ingly enough, some exact solitonic solutions were recently 
found for this equation in one dimension, but with time- 
and space-dependent interatomic interactions and trap- 
ping potential strengths 13, U, 13, 13, 14 1. Obtaining 
such exact solutions allows for testing the validity of the 
Gross-Pitaevskii equation at high densities, obtain the 
long-time evolution of the soliton where numerical tech- 
niques may fail, and helps to understand soliton forma- 
tion and propagation. 

The so-called Darboux transformation method [3] was 
used to obtain such exact solutions. We have shown in 
a previous work [191 ] that such exact solutions may exist 
only for specific functional forms for the interatomic in- 
teraction and trapping potential strengths. For example, 
in the work of Liang et al. [1 21 ] , the exact solution is found 
only when the trapping potential is quadratic, expulsive, 
and the interatomic interaction strength is growing ex- 
ponentially with time with a rate that equals the trap- 
ping potential strength itself. Such a restriction, makes 
the exact solution less interesting from an experimen- 
tal point of view. In an attempt to soften this restric- 
tion, we found that exact solutions may also be found 
for constant, linear, or quadratic potentials, and with 
interatomic interaction stren gth s that can be constant, 
growing, or decaying in time [191 ] . 

Here, we exploit our previous result to obtain ex- 
act solitonic solutions of the time-dependent Gross- 
Pitaevskii equation with linear trapping potential and 
constant interatomic interaction strength that can be 
positive or negative. The Gross-Pitaevskii equation de- 



scribes, in this case, the surface of the condensate. 

The Gross-Pitaevskii equation. — Near the surface of a 
Bose-Einstein condensate, the quadratic trapping poten- 
tial can be approximated by a linear potential and the 
surface of the condensate can be regarded as an infinite 
plane. The Gross-Pitaevskii equation, in this case, takes 
the form Ei SI 
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where ip(x,t) is the condensate wavefunction, and x is 
the coordinate normal to the surface of the condensate 
such that the bulk of the condensate exists in the region 
x < 0. The force constant F arises from linearizing the 
harmonic potential near the surface, namely F — mu) 2 R, 
where u> is the characteristic frequency of a spherically 
symmetric harmonic trapping potential, R is the radius 
of the condensate, and m is the mass of an atom. The 
interatomic interaction strength is proportional to the 
scattering length a which can be positive or negative. 

The characteristic length 5 in the surface region is de- 
fined by equating the kinetic energy fi 2 /2m5 2 to the po- 
tential energy FS, namely 5 = (h 2 /2mF) 1 / 3 . Scaling 
length to S, time to r = 2mS 2 /%, and the wavefunction 
to y/po = l/^8n\a\S 2 , the previous equation takes the 
dimcnsionlcss form 
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where we have retained the nonscalcd symbols for conve- 
nience. The parameter p 2 = — sgn(a) allows for treating 
the repulsive case (p 2 = — 1) and attractive case (p 2 = 1) 
simultaneously. 

For the case of repulsive interactions, the Thomas- 
Fermi approximation can be used to estimate S and po 
in terms of the size of the condensate, R, and the central 
density ptf [13]. It turns out that S/R w 7~ 4 / 15 and 
Po/Ptf ~ 7~ 4 ^ 15 , where 7 — Na/ao is the dimensionless 
interaction strength, ao = y/Ti/mu) is the characteristic 
length of the harmonic oscillator potential, and is the 
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number of atoms. For a typical Rb condensate with 
10 4 atoms the two ratios are roughly equal to 1/4. The 
unit of time is, in this case, given by r w 27~ 1 / 15 /ui. Ap- 
proximating the quadratic potential of the Bose-Einstein 
condenste by a linear one in the surface region is accurate 
only within a region of width S around the Thomas-Fermi 
surface. It turns out, however, that some of the solitonic 
solutions we obtain here have width larger than ft Fur- 
thermore, the dynamics of these solitons is such that they 
drift from the surface region towards the bulk region of 
the condensate where the potential is not linear anymore. 
Therefore, such solutions can be considered only as the 
initial states of the time-dependent solitonic excitations 
of the condensate. 

The Darboux Transformation and the New 
Solutions. — The first step in the Darboux trans- 
formation method is to find a linear system of equations 
for an auxiliary field *&(x,t) such that Eq. ([2]) is its 
consistency condition [3] • Using the method described 
m Ref. 191, we find that the following linear system 
corresponds to Eq. ([2|) 

*, =J*A + U*, (3) 
i^t = W* + 2(CJ + U)*A + 2J^A 2 , (4) 

— •<*.<> - ati itl> - 
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I 
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W = (C 2 - z/2)J + 2CU - J(U 2 - U x ), C(i) = it/2, 
and Ai and A2 are arbitrary constants. The subscripts 
x and t denote partial derivatives with respect to x 
and t, respectively. Equation ^ is obtained from the 
consistency condition ^ xt — ^tx and by substituting 
q(x, t) — r*(x, t) = ip(x, t). 

This linear system of 8 equations, Eqs. ((3]) and f4|, 
reduces to an equivalent system of 4 equations with non- 
trivial solutions by making the following substitutions: 
Ai = — A|, <pi = -02j ^1 = -P 2( f>2- Tne reduced 
system can be solved once the so-called seed solution 
r(x,t) = q*(x,t) is specified. We have also shown in 
Ref. [19( that the wavefunction ip (x,t) = Aexp(£<^>o) 
with (po = t(p 2 A 2 — (i 2 /3 + x)), where A is a real con- 
stant, is an exact solution of Eq. ([2]), and thus can be 
taken as the seed solution. 

The Darboux transformation can now be applied to 
the linear system to generate a new solution of Eq. @ 
as follows [Si 

^(x,t)=M^t)-^(X 1 +Xl)cf >2 r2/(p 2 \<t>\ 2 + m 2 )- (5) 

Substituting for ipo(x,t), ip2{x,t), and 4>%{x, t), we obtain 
the following new exact solutions to Eq. |2]): For the 
repulsive interactions case {p — ±i), the solution is 
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2u+ cosh 6»- 2iw+sinh6> + (|u+| 2 + 1) cos (3 + i(\u+\ 2 - l)sin/3 



(|?i+| 2 - l)sinh6» + 2u+sin/3 
and for the attractive interactions case (p = 1), the solution is 

2u+ cosh6>- 2iuf sinh6» + (|u+| 2 + 1) cos/3 + z(|u+| 2 - l)sin/3 



i/j(x,t) 



A - V8X lr 



(\u+\ 2 + l)cosh0 + 2i^ cos/3 
I 



(6) 



(7) 



where 9 = \[2 [A r (t 2 + x) + 2(A r X u - A { X lr )t] - S r , 
f3 = -V2 [A,(t 2 + x) + 2(AjAii + A r X lr )t] + ft , tt ± = 

V8 P A/b ± , b ± = 4X1 ± A, A = ^2Xf-p 2 A 2 , and 5 
is an arbitrary constant. Here, the subscripts r and i 
denote real and imaginary parts, respectively. 

Properties of the Solutions. — Here, we describe the 
main properties and features of the exact solitonic solu- 
tions found above. The five arbitrary constants 6 r , ft, 
Aii, Ai r and A, have the following effects on the solu- 
tions: The constants ft and ft have the trivial effect of 
shifting the solutions in the x- and ^-coordinates. There- 
fore, we set from now on ft = ft — 0. The other three 



constants are combined under a square root in the above 
expression for A. Thus, depending on the values of these 
constants, A can be real, imaginary, or complex. If A 
is real, (3 will have no ^-dependence and the solutions 
will be nonoscillatory, i.e., single-soliton solution. If A 
is imaginary, 9 will have no x-dependence, and the so- 
lution is oscillatory. If A is complex, the solution will 
be a combination of both previous cases, namely oscilla- 
tory but with a localized envelope which has a density 
profile that is similar to that of gap solitons (See Fig.l 
of Ref. These cases are shown in In Table 1, and 

specific cases are visualized in Fig. [TJ where we plot the 
density p(x,t) — |^(a;,t)| 2 at t = for the attractive 
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FIG. 1: Density p(x) = |V>(a;)| 2 at time t = for the case 
of attractive interactions. The arbitrary constants chosen to 
generate these plots are: S r = Si — and A = 1 for all plots. 
In (a) Aii = 0, Ai r = 0.29, in (b): Xu = 0, A ir = -0.6, in 
(c): Aii = 0, Ai r = 0.8, in (d): Ai, = 0, A ir = —1.5, in (e): 
Aii = 2, Air = -0.6, and in (f): Xu = 2, Air = 0.29. Typical 
values of the length and density units, 5 and po, are given in 
the text. 





Air > 


Air < 


A u = 


multi-solitonic 
with broad edges 


multi-solitonic 
with sharp edges 


Aii < 


single-solitonic 
with broad edges 


single-solitonic 
with sharp edges 


X u > o 


multi-solitonic 
with envelope 


multi-solitonic 
with envelope 



TABLE I: Classification of the solitonic solutions. Fig. la 
shows an example of the multi-solitonic solutions with broad 
edges, Fig. lb shows an example of the multi-solitonic solu- 
tions with sharp edges, Fig.lc shows an example of the single- 
solitonic solutions with broad edges, Fig. Id shows an example 
of the single-solitonic solutions with sharp edges, and Fig. Id 
and le show an example of the multi-solitonic solutions with 
an envelope. 



interactions case. 

For the special choice Xu = oo, the coefficient u + = 0, 
and the solutions , Eqs. (O and (?}, reduce to the simple 
forms: 



ip(x, t) = e**° (A ± iV8X lr e^cschj) 



FIG. 2: Color online Surface plots of the density p(x,i) = 
\i/>(x, t)\ 2 versus x and t for the attractive interactions. The 
upper plot corresponds to Fig. 1(d) while the lower figure 
corresponds to Fig 1(a). 
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These two equations show that our solution consists of a 
bright soliton embedded in the background 

The trajectory of a given soliton peak is obtained 
from the condition 9 — (for single solitons) or (3 — 
(for multiple solitons). The former gives x = —t 2 — 
2(Aij — Ai r Aj/A r )< + S r /y2 and the latter gives x = 
-t 2 - 2{X U + X lr A r /Ai)t + 5i/V2. The trajectory is 
thus parabolic in time with an acceleration of -1. In real 
units, this acceleration equals —F/m which is equal to 
the acceleration associated with the gravity-like waves 
propagating on the surface of the condensate [2^ 
This behavior can be clearly seen in Fig. [2] 

Using the above relation between x and t, we can elim- 
inate x from p to obtain the peak soliton density as a 
function of time, i.e., along the trajectory. This is shown 
in Fig. [3] where we notice that the soliton peak oscillates 
between a minimum and a maximum. The frequency 
of the oscillation equals v8| A| 2 Ai r /Ar. For the case of 
Fig- [He), where u+ > 0, the maximum appears at times 
defined by (3 — (2n + l)ir, namely t = mrA r /V8\ A| 2 Ai r , 
and the minimum appears at times defined by (3 — 2mr, 
n = 0, 1, 2, .... The maximum peak density is given 
by (A + \ASAir) 2 and the minimum peak density is given 
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FIG. 3: The soliton peak density along its trajectory for 
attractive interactions. The upper figure corresponds to 
Fig. 1(c) and the lower figure corresponds to Fig. 1(d). 



by (A - V8Air) 2 . For the case of Fig. QJd), the situ- 
ation is reversed since < 0. During this peak os- 
cillation, the number of atoms in the solitons is being 
exchanged with the background maintaining a dynamic 
stability 12j. The frequency of atoms exchange is con- 
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stant with time, unlike the case of Ref. 

An interesting feature of the soliton is its phase. We 
plot in Fig. [4] the phase of a multiple-soliton solution. 
This shows that the phase difference between the main 
neighbouring solitons is ir. It shows also that this phase 
difference, the width of solitons, and the distances be- 
tween them do not change with time. The soliton trains 
realized experimentally in a one-dimensional condensate 
with attractive interactions Q are very similar to our soli- 
tonic solutions in Fig. 1(d) and (e). It should be noted 
here that this is the case since, in addition to describ- 
ing the surface of the condensate, Eq. O describes also a 
one-dimensional condensate as long as the density is not 
too large [25} . The center of mass motion of the soliton 
train in our case is different than in the experiment of [5] 
due to the fact that we use a linear potential while in the 
experiment a quadratic potential is used. On the other 
hand, the density profile, relative phase, and number of 
solitons in the train can, in principle, be accounted for. 
This is the case since, in our theory and in the experi- 
ment, the width of the individual solitons (of order 5) is 



FIG. 4: The density (solid curve) and phase (light curve) of a 
solitonic solution for attractive interactions with S r — Si — 0, 
A — 1, Xu = 0.03, and Ai r = 0.6. The upper plot is for t = 
and the lower plot is for t = 0.85r. 



much less than the that of the background (of order R) . 

Since the density profile of the exact solitonc solu- 
tions found here depend only on the coordinate per- 
pendicular to the surface of the condensate, we predict 
that in a spherical condensate, a three-dimensional shell- 
like soliton may exist which in a sense similar to the 



ones reported in Refs. [26|,|28|]. Furthermore, the above- 
described dynamics indicates that, starting from the sur- 
face, this shell will be shrinking in radius. 

For the case of repulsive interactions, Eq.[S] shows that 
the density diverges at certain points along the parabolic 
trajectory described above. These points correspond to 
the points in Fig[3]where the soliton peak density is max- 
imum. The fact that the density diverges, does not make 
this solution nonphysical, since the number of atoms in 
the soliton is finite. 

It should be mentioned that exact solitonic solutions of 
Eq. @ have been essentially obtained using the so-called 
inverse-scattering method [27] . However, the present 
work represents another method of obtaining such exact 
solitonic solutions. We believe that our method in [lj| of 
obtaining the Lax pair is more systematic since the Lax 
pair of Ref. [27j was introduced as an assumption. Fur- 
thermore, in Ref. 27j, only formal solutions are derived 
for the soliton train case. This is in contrast with the 
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present work where we obtain explicit single as well as 
multiple solitons solutions. Finally, while in Ref. 27j . 



only the attractive interactions case is considered, we 
have derived here solutions for both the attractive and 
repulsives cases. 

In conclusion, we have found exact solitonic solutions 
of a time-dependent Gross-Pitaevskii equation with lin- 
ear trapping potential for both cases of attractive and 
repulsive interatomic interactions. These solutions may 
be regarded as solitons in the surface region of a three- 
dimensional Bose-Einstein condensate, or solitons in a 
one-dimensional condensate. 
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